Complexity perspectives: an anomalous diffusion approach 



OS 

o 
o 

(N 

as: 
< 



Guilherme R. Rezende 1,2 Q Luciano C. Lapas 1,3 |l| and Fernando A. Oliveira^ 

1 Instituto de Fisica and Centro International de Fisica da Materia Condensada, 
Universidade de Brasilia, Caixa Postal 04513, 70919-970 Brasilia, Distrito Federal, Brazil. 
2 Escola Tecnica Federal de Brasilia, 70730-521 Brasilia, 
Distrito Federal, Brazil. 3 Departamento de Matemdtica, 
Universidade de Brasilia, 70910-900 Brasilia, Distrito Federal, Brazil. 

The science of complexity is far from being fully understood and even its foundations are not well 
established. On the other hand, during the last decade, the random motion of particles or waves 
- the so-called diffusion - has been known better. In this paper, we discuss how simple ideas of 
diffusion can be used to deal with the description of most complex structure. 
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ANOMALOUS DIFFUSION AND COMPLEXITY 
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Complexity is still in its begin due to the large class 
of proposal and problems it reaches. Some few attempts 
have been done in order to measure complexity 0, [3] • 
However, complexity is "build up" from disciplines such 
as non-equilibrium statistical mechanics, non-linear dy- 
namics, and sciences as diverse as geophysics and biol- 
ogy, which are "complex" themselves, and consequently 
any progress in its branches will contribute to its con- 
ceptual evolution. Non-equilibrium statistical mechanics 
is somehow the simplest part of this high block, and its 
evolution shines light on all those connected area. Most 
developed part of the non-equilibrium statistical physics 
is concerned with relaxation process, i.e. the way a given 
system relax to equilibrium or to a steady state. Few the- 
orems have been established for that, mainly from linear 
response theory were one can compute the ansatz of a 
system to a small perturbation 

Complex systems are usually nonlinear, inhomoge- 
neous, and far from equilibrium. A common feature of 
those systems are pattern formations Q , where the trans- 
lation symmetry is broken. Pattern formations are com- 
mon as well in diffusion of alive beings [5|, |6| , and it has 
been demonstrated @, S] that the violation of ergodicity 
presents a situation similar to that of anomalous diffu- 
sion 0, ■ This is a green signal to extend those ideas 
to more complex systems. However, we need be careful 
to not drive ourselves in wonderful, but unrealistic gen- 
eralizations. Ergodic hypothesis (EH) is a fundamental 
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issue in statistical physics [llj, [12|, |13|, [14| and it validates 
a large number of theorems in physics. As a consequence 
it becomes a branch of mathematical physics, and a large 
number of works dealing with EH have been done by 
mathematicians, with small value for those interested in 
the analysis of mensurable quantities. In this context the 
Khinchin theorem (KT) is of special importance, since it 
was formulated in terms of correlations functions which 
can be obtained by use of linear response theory, and 
they are also directly connected with experiments, such 
as light or neutron scattering [3]. 



Diffusion is of fundamental importance for physics and 
much more simple to describe than reaction rates 0,0, 



181 ] , besides being the focus of extensive research in many 
interdisciplinary sciences d El M, ED, El HI E2, [Hf . 
An usual manner to study the diffusive dynamics is to 
deal with the mean square displacement of the particles, 
given by 



lim (x 2 (t)) cx t a , 



(1) 



where (. . .) denotes an ensemble average. The exponent 
a classifies the type of diffusion: for a = 1, we have nor- 
mal diffusion; for < a < 1, subdiffusion; and a > 1, 
superdiffusion. There is another way to deal with dif- 
fusion by considering a Brownian Motion described by 
a generalized Langevin equation (GLE) , or Mori formal- 
ism [2611 , of the form 



dP (t) 
dt 



U(t- t') P (<') dt' + F (t) 



(2) 



where P is the particle linear momentum, n (t) is the 
memory function, and F (t) is a random force, which ful- 
fills (F (£]} = and the fluctuation-dissipation theorem 
(FDT) [23: 



(F(t)F(t')) 







(3) 



In most of the experimental situations which the EH docs 
not holds arises in complex nonlinear or far from equi- 
librium structures where the detailed balance does not 
happens. Consequentl y, t he FDT does not hold too. 
In supercooled liquids [281 ] and systems with activeted 
dynamics (29| we can find good examples. Those sys- 
tems do not have apparently any easy analytical solution. 
Systems displaying glassy behavior have a multitude of 
metastable states that they visit after overcoming an ir- 
regular distribution of energy barriers, they exhibit aging 
and memory effects inherent to an activated dynamics 
characterized by a hierarchy of relaxation times [29( . On 
the another hand, diffusion can present closed solution 
for the main expectation values and arises as a simple 
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laboratory for discuss those properties. In this case, as 
we shall see, Lapas et al. give a full description of the 
KT validity [3(3]. In principle, it is generally possible 
to derive a GLE for Markovian systems by eliminating 
variables, whose effects are incor por ated in the memory 
kernel and in the colored noise [31[. Altogether, some 
results obtained for the GLE formalism should be valid 
for diffusion described by fractional Fokker-Planck equa- 



tions, since both formalisms yield similar results |32|,|3 



Khinchin theorem 

Since at the Khinchin times, most of the process stud- 
ied were driven to the normal diffusion with exponential 
relaxation. Thus, it is quite natural to inquire if the KT 
will be affected, for instance in the slow relaxation dy- 
namics that occurs in anomalous diffusion, where power 
laws, stretched exponential, Bessel functions and even a 
large jungle of functional behavior is possible for relax- 
ation, besides the velocities distributions may not be a 
Gaussian. 

Let us consider the Relaxation function 



R(t) = 



Cpjt) 

c P (py 



(4) 



where Cp(t) is the correlation function for the dynamical 
operator P, which is defined as 



C P (t)=(P(t)P(0))-(P(t))(P(Q)) 



(5) 



Explicitly the KT states that a system shall be ergodic 
in P as long as 



R(t -> oo) = 0. 



(6) 



In other words, irreversibility is a necessary and sufficient 
condition for EH. For macroscopic systems with a large 
number of degrees of freedom, the effect of past values of 
the forces usually vanishes for a sufficiently large t and 
the aforementioned condition is quite reasonable. 

Providing the well established method of recurrence 
relation 34|, Lee has attempted to the fact that irre- 
versibility is a broader concept than ergodicity and, in 
consequence, the KT may not work for all systems [11| . 
Physically, one can present some few systems where KT 
does not work, for example in the ordered phase of a mag- 
netic system Eq. © is fulfilled, nevertheless the system 
never goes through over all states, since the symmetry 
breaking allows only specific states. In this sense, Lee's 
statement creates a new challenge namely "in what sys- 
tems the KT works?" We shall show that the KT works 
for all range of anomalous diffusion described by a GLE, 
even in the presence of long range memory. The time 
averages of correlation functions are crucial for elucidat- 
ing the properties of dynamical processes and play an 



extremely important role in the ergodic theory and, con- 
sequently, in physics. For diffusive systems governed by 
the GLE we shall see that the condition R(t — ► oo) = 
is sufficient for the time average to be equivalent to the 
ensemble average, i.e., for the system to be ergodic. 

For any initial distribution of values, P (0) and 
(F(t)(0)) = 0, it is possible to obtain the temporal evo- 
lution of the moments of P from GLE, Eq. ([2|), 



(P(t)) = (P(0))R(t), 



and 



(P 2 (t)) = (P 2 ) eq +R 2 (t) <P 2 (0))-<P 2 ) 



eq 



(7) 



(8) 



To study relaxation we need to know R(t), which can be 
calculated analytically in restricted cases, being obtained 
numerically most of the times. Consequently, one can de- 
scribe completely those average values by the knowledge 
of R(t). Higher order moments can be obtained, but is 
not our goal here. Since Eqs. ([7]) and ([8]) are sufficient 
to show the validity of the KT for diffusion. If condi- 
tion (J6j> , R(t — > oo) = 0, is valid, then the time evolution 
will produce the ensemble average with (P (t)) =0 and 
(P 2 (t)) = (P 2 ) eq , i-e. EH holds, then the KT is valid. 
Now we may ask in what situation the Eq. ([6]) is not 
valid. First, one should note that the long time behavior 
is associated with the small values of z in the Laplace 
transform. Indeed, from the final value theorem we have 



lim R (t) 



lim zR (z) . 

z^0 



(9) 



Here we shall nominate the Laplace transform of a generic 
G{t) by G(z). Consequently it is only necessary to know 
R(z). Now taking the equation of motion, Eq. ([2]), multi- 
pling by p(0) and taking the ensemble average, we obtain 



dR(t) 
dt 



n(t-t')R(t')dt', 



(10) 



which Laplace transform yields R(z) = l/[z + H(z)]. 
Morgado et al. [2(| have shown that if 

n(* ->■()) « cz v , (ii) 

where c is a positive constant, then a = 1 + v. Con- 
sequently, the behavior of the memory n (z) for small z 
determines the long range behavior of the diffusion, Eq. 
(J). Now the limit into Eq. © reads 



lim R (t) 



lim I 



l + cz c 



(12) 



For most of the diffusive process < a < 2, the irre- 
versibility condition R (t — ► oo) = occurs. However, 
this condition fails for ballistic motion, a — 2, when 
R(t — ► oo) = 1/(1 + c) and the time-correlation function 
will be non-null for long times. Thus, if the ballistic sys- 
tem is not initially equilibrated, then it will never reach 
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equilibrium and the final result of any measurement will 
depend on the initial conditions. In this situation the 
EH will not happens, however once again the KT works 
since the violation of the EH was induced by the viola- 
tion of the irreversible condition, Eq. ©, as predicted by 
Khinchin. In other words, the ballistic diffusion violates 
ergodicity and the FDT [9j. The major consequence of 
the violation of condition, Eq. (J6j> , is the presence of a 
residual current. Let us suppose that the system starts 
with an average current such that (P (0)) 7^ 0. From 
Eq. 0, we obtain that for R (t — > 00) 7^ a residual cur- 
rent remains. However, the effective current can be very 
small compared to (P (0)) and its value, as any other 
measurable property for ballistic diffusion, will depend 
on the value of c. In other words, the system decays to a 
metastable state and remains in it indefinitely, even in the 
absence of an external field. For a > 2, R(t — >oo) = l 
, no memory of the initial conditions is lost, the pro- 
cess is not a diffusive one, and it is an activated pro- 
cess for which the GLE does not work 0]. Therefore, 
the results of this letter apply to all kinds of diffusion, 
< a < 2, described by a GLE independent of the mem- 
ory range, which gives origin to new studies in ballistic 
diffusion pj| [11 [IB] . 

For stationary systems, x(t,t') — x(t ~ one can 
define the time average integral as 



I ta = lim ^ 

T->oo T 



JQ 



X{t,t')dt'dt, 



Integrating by parts, we arrive at [111 ll 

1 



lim 

t — >oo 



t x ft 

X(t')dt' + R(t) - - / R(t')dt' 

I! WO 



(13) 



(14) 



Given that R(t) is a real-valued function that converges 
asymptotically to a finite value, since we are dealing with 
the linear momentum autocorrelation, we can use a gen- 
eralization of the final-value theorem for Laplace trans- 
forms 135], 



the irreversibility condition 
end up with 



holds. From Eq. ([9]) we 



X(0) = Xs 



(17) 



Again this is a consequence of the irreversibility condi- 
tion. Therefore, irreversibility is a necessary and suffi- 
cient condition for the EH to hold in diffusive processes 
described by a GLE. 

For nonlinear Hamiltonian maps [36| there is no gen- 
eral framework to address the problem and, in particular, 
the absence of a coupling to a thermal bath (explicit in 
the GLE) and consequently the lack of a detailed bal- 
ance relation or FDT may require a specific analysis of 
each case. However, since it is possible to give a kinetic 
description of the Hamiltonian dynamics by means of a 
fractional Fokker-Planck-Kolmogorov equation 37| , it is 
expected that the treatment of anomalous diffusion in 
such systems should also be possible by the GLE formal- 
ism. Further research in this direction is needed and will 
open new perspectives. 



Nonexponential behavior 

From the above results, it is quite clear that the corre- 
lation for ballistic diffusion will not decay exponentially 
to equilibrium. Besides the ballistic case, any anoma- 
lous regime will present nonexponential decay. Even for 
normal diffusion, a large number of relaxations may be 
nonexponential [20(. There are a large number of phe- 
nomena where the systems do not relax immediately to 
equilibrium. Those phenomena, usually associated with 
non-aging, have nonexponential relaxation and are most 
commonly described by power laws or stretched exponen- 
tials. The study of anomalous relaxation has produced 

" Mill El El, El. 



quite interesting results [28j, |3S 

For a system described by a GLE of the form Eq. ([2]), 
the evolution relies on the noise that drives the particles. 
For a harmonic noise |20| 



lim zR(z) 



lim ± f R(t)dt. 



From the above relation, we obtain from Eq. (|14p 

X(0) + R(t -> 00) - lim zR{z) = Xs, 
2^0 



(15) 



where Xs is the time independent value, the so-called 
static susceptibility. On the other hand, taking the 
Laplace transform of Eq. ([TO)) , we obtain x(z) + zR(z) = 
Xs- Taking the limit z — > 0, the previous relation be- 
comes 



X(0) + lim zR(z) 

z— *0 



Xs 



(16) 



Comparing Eq. (|T6|) with Eq. JT5J), one should conclude 
that the EH can only be valid if R(t — > 00) = 0, i.e. if 



F(t) = 



1 



V2fcBT 



\J p{u}) cos [uot + 4>{uj)] did, (18) 



where p is the noise spectral density, fc^ is the Boltzmann 
constant, and </> is a set of random phases in the range 
< 4> < 2tt. 



A systematic study carried on by Vainstein et al. [44 1 
has shown that the spectral density plays a fundamen- 
tal role in the description of stochastic processes as we 
shall see. First, the memory function H(t) can be easily 
obtained by the use of the FDT, Eq. ©, as 



n(*) 



p(u)) cos(ujt)dLu, 



(19) 



in such a way that the average cancels the random terms 
and obviously the memory is a deterministic even func- 
tion. The Laplace transform of the memory H(z) is an 
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odd function in z; therefore, R(z) is also an odd function 
in z. This implies that by inverting the Laplace trans- 
form, R(t) is an even function of t 



R(-t) = R(t). 



(20) 



Second, the condition given by Eq. to determine 
the exponent a (a = 1 + v) , for a spectral density of the 
form 



P (lu) = 



ato" , for lo < lu s 
<?(w), otherwise. 



(21) 



Here, the function g(u>) is arbitrary as long as it is suffi- 
ciently well-behaved and that its integral in the memory 
function converges. If one is interested only in the long 
time behavior t 1/lj s , it can be taken to be 0. With 
this noise density of states, it is possible to simulate many 
diffusive regimes 45[ . Noise of this form can be obtained 
either by formal methods or empirical data. Using this 
expression in Eq. (|19[) . taking the Laplace transform in 
the limit z — > 0, we have 



n{z) (x 



-az ln(z), 

z, 



for /3<1, 
for /3= 1, 
for (3 > 1. 



(22) 



Consequently, for this type of noise, there is a maximum 
value of a, i.e., a < 2 for any value of [3. It should be 
noted that the case (3=1 does not lead to a memory 
whose Laplace transform is in the form of Eq. (jlip . For 
— 1 < j3 < 1, we obtain a = 1 + f3. For j3 > 1, one has 
a = 2, which shows that ballistic diffusion is a limiting 
case for the GLE with this type of memory. In this way, 
we can emphasize that 



(3, (3<1; 
1, (3>1. 



(23) 



This shows that a noise spectral density in the form of 
a power law can produce only diffusive motion in the 
region < a < 2. Diffusive motion beyond ballistic is 
not allowed, what can be observed by using the Laplace 
transform. 

A recent study shows that exponential decay, power 
laws or even Mittag-Leffler functions are particular cases 
of a more general function [4J] which approximates the 



decay. Indeed, considering time reversal symmetry [15 1 
and the discussion above (see Eq. (|20[) L the correlation 
function must be even and cannot be any of those forms. 
We shall expose here the conditions under which it is pos- 
sible to obtain an approximately exponential decay, what 
happens in certain circumstances for normal diffusion. In 
this case, 



7 =limn(z) = -p(0). 
z— >o z 



(24) 



Consequently, the friction in the usual Langevin equation 
is nothing more than the noise spectral density for the 
lower modes. For an arbitrary memory, the system has a 
rich behavior; even for normal diffusion it is possible to 
show the existence of at least three time ranges 44|. A 
normal diffusion can be obtained using a spectral density 
of the form [H 



= J 7T 



27 



LO > W, 



0, 



(25) 



For a broad band noise spectral density r y/uJs < 1, and 
long times t > 7 -1 it is possible to decouple Eq. (fTU|) 
to obtain an exponential decay of the form R(t) — 
exp (— jt), which will bring the system to equilibrium, 
that is R(t — ► 00) — > 0. 



Evidencing the KT 

In order to illustrate the analytical results, we have 
numerically integrated the GLE, Eq. ([2|), to obtain ap- 
proximations to the probability distribution of particle 
velocities using histograms. 

In Fig. ([I]) we show the probability distribution func- 
tions obtained for subdiffusion ({3 — —0.5), normal diffu- 
sion (/? = 0), and superdiffusion ((3 = 0.5) for the values 
a = 0.25 and g(w) = 0. We have used ui s — 0.5 for 
all cases except for subdiffusion, which demands a broad 
noise lo s = 2 to reach the stationary state. In all cases, we 
expect that R(t — * 00) = 0, and that the EH will be valid 
even for the subdiffusive (superdiffusive) case, despite the 
fact that W = (W = 00). These relations can be seen 
by considering the limit W — lim z ^o R{z). If the EH 
is valid, the velocity probability distribution will be the 
same for an average over an ensemble of particles and for 
a time average over the trajectory of a single particle for 
long times after the system has reached an equilibrium 
state. Note that despite the presence of large fluctua- 
tions in the time average case due to numerical errors, 
there is a good agreement between the resulting ensem- 
ble and time distributions. The three probability distri- 
butions converge toward the Maxwell-Boltzmann distri- 
bution, which is in accordance with previous analytical 
results lldl. 



REMARKS 

In order to understand cellular life processes like 
cell division and nerve cell growth, it has been shown 
that the protein reaches the position via a random, 
one-dimensional diffusion movement along the micro- 
tubule [ifjj ]. Besides that, statistical analysis leads to 
purely diffusive dynamic laws such that some internal 
processes induce a breakdown of the central-limit theo- 
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FIG. 1: (Color online) Numerical results for the probabil- 
ity distribution function for subdiffusion (top, a — 0.5), nor- 
mal diffusion (middle, a — 1) and superdiffusion (bottom, 
a — 1.5). The time averages (circles) are obtained by follow- 
ing one particle trajectory and calculating the histogram for 
times from t = 100 to t = 5000. For the ensemble averages 
(squares), we calculate the histogram using 5- 10 4 particles, at 
time t = 1000. The continuous line is the Maxwell-Boltzmann 
distribution. Insets: Curves a correspond to the functions t a 
and curves b to the simulated mean square displacements. 



rem and leads to anomalous diffusion |47[ ■ Those molecu- 
lar motors are trapped in a rough energy landscape sub- 
jected to a kind of glassy behavior common to a great 



varietyof complex systems such as glasses [48[ and pro- 



teins 4t|, 5(J, just to mention some examples, which cor- 
responds to apply an external force strong enough to 
place the system near the stall point of the motor [5l| . 

A large number of theorems in Statistical Physics and 
even a proper applications of formalism such as a " simple 
linear response theory" relay under the ergodic hypoth- 
esis and its validation. The Khinchin theorem put the 
EH in a straight way, and it has a practical character, 
since it is exposed in therm of response function. In this 
way, the Lee's works need a large discussion concerning 
the validity of the KT in specific systems. In this work 



we have discussed the KT for anomalous diffusion, which 
are ergodic in the range of exponents < a < 2, where 
a defines the asymptotic behavior of the diffusion, Eq. 
((T|). For a — 2 we have the special ballistic case as we 
have seen the EH does not work for ballistic diffusion, 
however, the condition Eq. ^ does not work either, con- 
sequently the KT remains as a strong reference, at least 
for diffusion. In recent years Molecular motors has been 



driving a lot of attention [22J, |25| with a large potential 



for pure and applied science, there those discussion may 
be very useful. There are many situations of violation of 
the EH, particularly in glassy syste ms 13811 . and another 
were EH holds, for example recently [391 ] dynamical sim- 
ulations and equilibrium statistical mechanics were used 
to treat glass transition calculations. The agreement be- 
tween them is a strong indication of the validity of the 
HE. This is a quite surprising result for such complex 
system. Disordered systems are a large universe to ex- 
plore the basic assumptions of statistical mechanics in 
particular of the KT. 

This work was supported by FAPDF, CAPES, and 
CNPq. 
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